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Abstract 

In this work are discussed possible generalizations of the expanding locally 
anisotropic metric ansatz with respect to approximately Newtonian many body 
gravitational systems. This ansatz describes local point-like matter distribu- 
tions on the expanding Universe also allowing for a covariant parameterization 
of gravitational interactions at intermediate length scales. As an example of ap- 
plicability it is modelled a disk galaxy model matching the physical parameters 
of the galaxy UGC2885 and it is shown that, by fine-tuning the metric functional 
parameter, the flattening of the galaxy rotation curve is fully parameterized by 
this metric. It is further analysed the mass-energy density corrections due to the 
expanding anisotropic background being shown that although there are negative 
contributions within the galaxy plane the total mass-energy density is strictly 
positive. Outside the galaxy plane, by considering an anisotropic functional 
parameter, it is shown that the mass-energy density is also strictly positive. 



1 Introduction 



In this work our main objective is to generalize to many body gravitational 
systems the recently suggested expanding locally anisotropic metric ansatz [Tl[2] 
describing point-like local matter distributions in an expanding background, 
the Universe P], as well as to shown that such ansatz allows to parameterize 
corrections to the gravitational acceleration within intermediate gravitational 
scales. In particular, as an example, it is shown that for a specific fine-tuning 
of the functional parameter of this metric the flattening of galaxies rotation 
curves with respect to the predicted profiles, when considering only the Newto- 
nian acceleration due to baryonic matter within galaxies [UlS] , is fully described 
such that this gravitational background ansatz can be interpreted as a covari- 
ant parameterization of the gravitational corrections attributed to cold dark 
matter [6] . Assuming a thin exponential disk approximation we apply such pa- 
rameterization to a simplified lattice model roughly matching the characteristics 
of the galaxy UGC2885 [ZHlO]. 

The ELA metric ansatz was derived in [T] as a generalization of previous 
metric ansatze [TT] describing local matter distributions in the expanding Uni- 
verse [3] and noting that, although global space-time isotropy is mandatory, 
local matter distributions generate local anisotropy [HI [13]. Specifically the 
infinitesimal line element for this ansatz is 
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where U = 2GM/{c^ri) is the Schwarzschild (SC) gravitational potential, H = 
d/a is the time dependent Hubble rate defined as the rate of variation of the Uni- 
verse scale factor a, c is the speed of light in vacuum, M is the SC gravitational 
mass, G is the Gravitational constant, df2^ — rf{d6^ + s'm^ 9 dip) is the solid 
angle line-element and we are employing spherical coordinates {t,ri,ip,9) for 
which the radial coordinate is an area radius (the area of a sphere is A = 47rr^) 
such that the integration measure is independent of the Universe scale fac- 
tor, ^/—g = rfsm{6). The ELA metric ([1]) interpolates between the SC met- 
ric [13] for small radii values (near the SC event horizon at the SC radius 
Ti ~ ri se = IGM/c^) and the Robertson- Walker metric [TS] for large radii 
values (ri — >■ oo). The shift function depends on a functional parameter a 
which, at the SC event horizon must be greater or equal to a{ri,sc) > 3 to pre- 
vent this horizon from being a space-time singularity, at the origin must have 
at least a leading divergence by a{ri — >■ 0) ^ — 1/ri to ensure that the SC mass 
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pole value coincides with the SC mass M and, at spatial infinity, it must be 
finite such that the shift function asymptotically vanishes and the RW metric 
is recovered [T]. 

As for the mass of the extended gravitational background encoded in the 
ELA metric in excess to the mass of the cosmological background described by 
the RW metric it was shown in [TH] to be 



In this expression is the (extended) mass-energy contribution due to the 
expanding locally anisotropic metric in excess of the SC mass pole of value M 
and pbw — 3iJ^/(87rG) is the cosmological background mass-energy density. 
The value for this mass is finite either when it is considered a radial cut-off for 
a above which this functional parameter is null such that the isotropic McVittie 
metric II Ij is exactly recovered, or when the functional parameter is asymptot- 
ically null at spatial infinity being asymptotically proportional to ^ ^ with 
n> 2. Hence either of these limits for a must be considered to ensure a finite 
mass contribution due to the expanding anisotropic background. 

In addition, in between the two well established asymptotic limits (the SC 
metric and RW metric), there are no specific bounds on the functional parameter 
a such that it allows for a covariant parameterization of unniodelled observable 
effects at intermediate length scales. Also we recall that so far no direct phys- 
ical interpretation for the functional parameter a exist, here the ELA metric 
ansatz is interpreted at most as a covariant parameterization of gravitational 
interactions corrections compatible with the well established asymptotic back- 
ground solutions, the local SC metric and the RW metric. In this work we 
explore the application of such parameterization to describe the dynamics of 
many body gravitational systems, namely the fiattening of rotation curves for 
galaxies, which indicates that such parameterization may be physically inter- 
preted as cold dark matter. Such parameterization has also been applied to the 
Solar System dynamics |16j , in particular allowing for a parameterization of the 
variation of the Astronomical Unit [17] . 

2 Ansatze for many body systems 

With the objective of generalizing the ansatz ([1]) to many body gravitational 
systems let us note that, generally, the above ansatz for the ELA metric ([1]) 
is obtained from any given metric g^^i, describing the gravitational background 
of local matter distributions by considering the radial shift function N'^j^j^ = 
—Hn {^/gm)°'~^^ / c. Hence to generalize the ELA metric ansatz to many body 
backgrounds it is first require to define the many body SC metric background. 
However no explicit analytical solution for a metric describing such background 
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is known. The standard approach is to consider a perturbative post-Newtonian 
solution to a specific order in the gravitational field U [1^ . Here we will consider 
the Newtonian approximation which is commonly employed in macroscopical 
galaxy models [TH] (see also [IS]), hence we define the local Newtonian grav- 
itational background for N massive bodies to first order in the gravitational 
field 

dssc.N = (1 - Un) {cdtf - dr^ - dn^ , 
^ 2GM„ 

Un = / Un , Un — -IT. r , 

c2|ri-ri.„| 

where Un is the total SC gravitational potential being C/„ the SC gravitational 
potential for each of the many bodies. The distance to each body position is 
given as usual by the Euclidean distance between the point at which the metric 
is being evaluated (ri) and the position of each massive point-like body (ri.n)- 
We recall that for simulations within the Solar System the gravitational field is 
at most of order U < 10~* and a second order expansion on the gravitational 
field is commonly employed, hence being obtained an accuracy on the metric 
definition of order ~ 10~^^ which is enough to match the current experimental 
accuracy for observational data }21j . In this work we will consider a galaxy model 
based in an analytical thin disk approximation for which the SC gravitational 
field is typically of order J7„ < 10~^ such that the second order corrections 
would be at most of order ^ 10^^^. Hence the first order expansion considered 
for goo.N and grr.N is accurate to order ^ 10^'' being enough for exemplification 
purposes and to obtain a macroscopical estimative for the functional parameter 
of the EL A metric ansatz. A more detailed model would require a higher order 
expansion on the gravitational field. 

Given the many body perturbative background ^ it is straight forward 
to generalize the ELA metric ansatz by considering the definition of the shift 
function with respect to the total gravitational potential 

ds'ELA.oS = {l-UN)icdt)^-dn^ (4) 

- {dn + N^ELA.efficdtyf . 

This metric has the same properties of the ELA metric ([l} , in particular at the 
event horizon for each of the n bodies there is no space-time singularity. We re- 
mark that the functional parameter a^ff parameterizes a collective gravitational 
correction to the SC metric due to the many bodies in the gravitational system 
such that a direct fit to the galaxy rotation curve can be carried allowing to 
map this functional parameter within the galaxy plane, hence to compute the 
mass-energy density and anisotropic pressures within the galaxy. 

To further compute the mass-energy density and anisotropic pressures out- 
side the galaxy plane, it is required to estimate the functional parameter acff 
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outside the galaxy plane. To achieve such estimate it must be defined each in- 
dividual body contribution in the gravitational system to the metric shift func- 
tion. Also aiming both at describing gravitational systems by more detailed 
models including the individual known massive bodies |17| , as well as to achieve 
a physical interpretation for the functional parameter of the ELA metric, it is 
desirable to have a functional parameter for each individual body in any given 
gravitational system. Hence, next, we are mapping the shift function N^j^j^ 
into an equivalent product of the many body contributions defined with respect 
to the gravitational fields of each individual body in the gravitational system. 
To achieve such factorization let us note that the event horizons for each of 
the many body is no longer a spherical surface, instead the Schwarzschild sur- 
face for each of these n bodies is defined by the 2-dimensional solution to the 
equation Um = 1 in the neighbourhood of each of the point-like massive bodies. 
Specifically, for a given body n, the event horizon is the solution of the following 
equation 



The solution of the radial coordinate ri,sc.n for each of the N event horizons 
can be parameterized by the angular coordinates 9 and Lp such that, for each of 
the n bodies, this equation can be expressed in terms of a modified gravitational 
potential as [/„ = !. Specifically, this modified gravitational potential, is defined 
as 



For gravitational systems for which the values of the Schwarzschild radius is 
much smaller than the distance between bodies, |ri.n — ri.gc.n(^, ^ I^Li ^ rj.nC^, "ys)!, 
this modified gravitational potential approximately matches the SC gravita- 
tional potential [/„ ~ C/„. 

Once we have defined these modified potentials, the shift function for the 
ELA metric can be factorized into a product of exponentials (1 — Un) which 
exactly vanish at the event horizon corresponding to each of the many body 
in the system and asymptotically converge to unity at spatial infinity, hence 
having the same properties of the original single body ELA metric. Specifically 
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we obtain the ansatz 



^ n=l 
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The map between the functional parameters of metrics (O and is straight 
forwardly obtained to be 



N 

(an + 1) hl(l - Un) 



With respect to the gravitational interactions corrections, these two parameters 
have a distinct meaning. While the functional parameter q;„ describes the grav- 
itational interactions due to each body n in the remaining iV — 1 bodies, the 
functional parameter deft describes the local effect near each body n due to the 
gravitational interactions of the remaining iV — 1 bodies. 

Most astrophysical systems, such as galaxies, lye approximately on a 2- 
dimensional plane. Hence the physical laws derived from direct observational 
data of a given astrophysical system are valid within such plane. The analysis of 
physical laws in the neighbourhood of the system is usually inferred from known 
gravitational laws as well as by analysing indirect observational data such as the 
gravitational lens effect. Following this discussion, from the map ([S]), known a 
profile for the functional parameter it is possible to derive a profile for 
each of the functional parameters a„ aiming at obtaining a estimation for the 
gravitational corrections due to the ELA metric outside of the plane of the 
astrophysical system. 

For a given gravitational system we are considering the following simplifica- 
tion assumptions: 

• the system is approximately planar allowing for a circular symmetric pla- 
nar model; 

• each of the functional parameters is circular symmetric with respect 
to the centre of mass ri.„ of each body M„ in the plane of the system; 

Given these assumptions the map ^ allows to obtain a scaling law for the 
functional parameter q;„ with respect to each body mass M„. Let us consider 
Nm bodies, each with a mass Mq, hence with total mass M = Nm Mq. If 
we consider these Nm bodies to be infinitesimally close to each other at a 
distance ri.o from a given test mass such that f7i„ = JJq ^ 2GA/o/(c^ f'l.o) 
and Uij^ — NUq (for i^r = 1, . . . , N), the effective functional parameter ajy^^j is. 
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according to map ([5]), unj^^ = — l+A^(ao + l)- Further noting that Nm = M/Mq 
and that the Nm bodies must be consistently described by one single body of 
mass M we obtain the scaling law for the functional parameters a„ with respect 
to some reference mass Mq to be 

M 

an = -l + ^K + l). (9) 

In particular, assuming that this scaling law is circular symmetric, allows to 
obtain a reference profile for the functional parameter ao on the planar system 
from a known profile for UcS, also on the planar system. Further assuming 
either approximately spherical symmetry for each functional parameter an or 
by parameterizing its anisotropy along the orthogonal direction to the planar 
system the functional parameter acs can be computed outside the planar sys- 
tem. For a given many body system, we will explore this relation to estimate a 
3-D map for acg. We also remark that this scaling law implies relatively large 
values of the parameter ckj\/ for large masses as it scales linearly with the mass. 

In addition, assuming the above scaling relation for each of the many bodies 
with masses M„ , in the limit of large radii the map ([SJ is approximately given 

by 

aoff(''i +00) 



Mat 

where Mjv is the total mass of the gravitational system, the masses Mn{0,(p) 
are defined in ([SJ and the last expression for acff is obtained by considering 
the approximation M„ « M„. Noting that far from the astrophysical body 
the gravitational potential approximately equals the gravitational potential of a 
single point-like massive body with the total mass of the N-body astrophysical 
system Un ~ 2G / {P ri) such that the mass contribution due to the ex- 
tended anisotropic gravitational background is expressed by equation ^ with 
M = Mm and further recalling that, as discussed in the introduction, to en- 
sure a finite mass contribution the metric functional parameter must either be 
asymptotically null at spatial infinity or strictly null above a finite radial cut-off, 
we conclude that the above asymptotic expression for aeff(''i ~ +00) must be 
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null. Hence, allowing for a generic value for the reference mass Mq, we obtain 
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These limits constitute a consistency check for the map ([S]) and the definition of 
the reference profile ao, in particular allow to set the asymptotic limit of such 
profile for large radii. 

We note that when considering a detailed model containing the known mas- 
sive bodies within a given gravitational system the ansatz ([7]) describing each 
individual body contribution to the background is applicable requiring a full 
numerical simulation with N functional parameters. However for a simplified 
model considering a planar average mass surface density, the ansatz ^ de- 
scribing one single effective functional parameter on the galaxy plane reduces 
significantly the complexity of the model allowing for a simpler parameteriza- 
tion of gravitational corrections. Hence in the remaining of this work we will 
employ the effective metric ansatz ^ to develop a simplified disk lattice model 
for galaxies similar to cosmological lattice models [24] . For a given galaxy, once 
a discretized 2D profile for the functional parameter aeff is obtained by fitting 
the rotation curve on the plane of the galaxy, the discretized profile for the 
functional parameter of each body on the galaxy model, a„ ([7]), can be param- 
eterized by the scaling law for um ® such that solving the system of equations 
corresponding to map ([5]) we obtain the discretized profile for ag on the galaxy 
plane corresponding to the model considered. From this discretized profile it 
is straight forward to explicitly compute the discretized profile for each of the 
a„, also on the galaxy plane, from the scaling law (|9]). Further assuming either 
spherical symmetry (or a parameterization of the anisotropy) of the parameter 
an with respect to the centre of mass of each body in the model the, value for 
the functional parameter aoS outside of the galaxy plane is computed by em- 
ploying the map ([5]) , hence allowing for a 3-D analysis of derived gravitational 
quantities in the neighbourhood of the galaxy such as mass-energy, anisotropic 
pressures and equation of state corrections due to the ELA metric background. 

When computing the equations of motion for cither the metric or ([7]), 
there is one more problem to address, although for the SC metric we have 
considered the series expansion to first order in the gravitational field for the 
metric components gsc.oo = 1 ~ Un and gsc.rr = 1, we have not considered a 
series expansion for the shift function N^j^j^. Although for small values of the 
metric functional parameter |a| ^ 10 an expansion of the factor (1 — U)" to 
the same order of the metric component 500 maintains the accuracy of the ELA 
metric approximation, for higher values of the functional parameter a higher 
order expansion of the shift function on the gravitational field U is required to 
attain the same accuracy. We recall that a similar problem is verified in the 
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PN formalism as the component 1700 requires an higher order expansion than 
the remaining metric components |18| . Although an higher order expansion of 
the ELA metric shift function can be considered, there is no clear technical 
advantage as it will difficult the derivation of the equations of geodesic motion 
and other derived quantities. Instead, for technical simplification we will carry 
the following computations to first order on the gravitational field with respect 
to the local metric background components gsc.oo ~ ^ ~ Un and gsc.rr = 1 
while considering the exact expression for the shift function N"^^^. Also this 
construction can be justified by noting that the shift function is here a functional 
parameter, hence a generic function for which the specific dependence on the 
gravitational field is unknown such that the specific order of its series expansion 
in the gravitational field cannot be exactly determined. In addition we will 
consider a series expansion for the metric and remaining derived quantities of 
second order on the Hubble rate H , hence of second order on the shift function 

Given this modelling setup we proceed to compute the metric connections 
and obtain the dominant contributions to the radial acceleration of a test mass 
on the many body background described by the ansatz (|4]). Considering that 
the many body system lies approximately on a plane of constant 9 and for 
non-relativistic velocities ri <C c, </ji <C c and 9i ^ c, the radial acceleration 
approximated to first order on the gravitational field Un and second order on 
the Hubble rate H is 



n if 



sinO 



2 (12) 



= (1 - C/at)"-'- ((1 - Un) (2 + ri log(l - C/jv)4ff) 

-2(1 + q){l - Un)^-^ - ri(2 + aeff - Un) t/^) , 

where dotted quantities represent differentiation with respect to the time coor- 
dinate t and primed quantities differentiation with respect to the radial coordi- 
nate ri . The acceleration component Fn is the standard gravitational Newton 
acceleration, F^ is the standard centripetal acceleration and Fjj2 is the lower 
order correction in the Hubble rate H due to the expanding locally anisotropic 
background as described by the ELA metric 

As long as the parameter aes is finite, for small values of ri, F^^ « is 
negligible, while for large enough values of ri the acceleration component F^^ 
is positive being approximated by the usual RW acceleration lim, j_>oo ^if^ « 
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Frw = —qH^ri, hence coinciding with the cosmological acceleration outwards 
the observer. In between these two asymptotic hmits, depending on the value 
of the functional parameter aes, the component Fj^i may either be positive or 
negative. We note that, by numerical inspection of , for large negative values 
of the parameter aes the correction to the Newton gravitational acceleration is 
negative, hence towards the central mass (TB] . In particular the metric functional 
parameter allows for a covariant parameterization of the excess acceleration 
towards the core of galaxies, hence increasing the orbital speed of massive bodies 
on the galaxy and allowing to describe the experimentally observed flattening of 
galaxy rotation curves |31[7]- In the next section we explicitly derive a discrete 
disk galaxy model approximately matching the physical parameters of the galaxy 
UGC2885. 

As for the mass-energy density for the many body ELA metric (|4]) it is 

= -^^^{l-UN)lN^ELA.cs{NELA.cS + 2ridrN^ELA.ce) 
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(1 + Un){1 - UNr^'-' i (1 - Un^x (13) 



X (l + ^log(l - UN)dra,s) - ^^((1 - a,s){deUNf 

1 - "off , 



■{d^UNf + 16nil-UN)drUN) 



In the limit of large radii, as long as acs is finite, this expression coincides with 
the RW mass-energy density p^w = 3H^/{8ttG). 



3 A Lattice model for galaxies 

In this section we will apply the many body ELA metric (U) to develop a simpli- 
fied lattice disk model for galaxies showing that the flattening of galaxy rotation 
curves can be fully parameterized by this ansatz. To model the galaxy bary- 
onic matter we are considering a discrete polar lattice on the galaxy plane with 
point-like masses at each lattice vertex. Although it is possible to implement 
an analytical model on the plane of the galaxy, to compute the values of the 
effective functional parameter acS ® outside the galaxy plane employing the 
method described in the previous section it is required to consider a gravita- 
tional system of point like masses as it employs a reference profile ao for each 
massive body (|9]). We also note that an explicit discretization of the mass-energy 
density simplifies the computational simulation of the galaxy model as we will 
develop in the remaining of this section. To explicitly compute the values of 
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the point-like masses on galactic plane we consider the same model employed 
in [5] consisting of an infinitesimal 2D thin exponential disk [19] describing the 
intergalactic gas and a 3D thin exponential disk [22] describing stellar matter 
across the galaxy plane plus a bulge matter distribution near the centre of the 
galaxy describing the galaxy core. The galaxy disk is considered to be finite 
having an explicit cut-off describing the galaxy edge [23]. 

For exemplification purposes we are modelling the galaxy UGC2885. Ex- 
perimental photometric data and red-shift measurements for this galaxy are 
available in the works [H[7yi0|. We note that within the analysis carried in 
these references there are discrepancies in the derived physical parameters of 
the galaxy. These discrepancies are mainly due to the distinct estimates for 
the Heliocentric distance D, which is derived from the systemic velocity Vgys 
which, in turn, is computed from the average of the measured red-shift. While 
in [3J|7] it is considered the uncorrected systemic velocity Vsys = 5794fcms~^ 
corresponding to a Heliocentric distance of £> = 118 Mpc, in [SJIH] the systemic 
velocity is converted to the motion relative to the local group and the cosmic 
microwave background Vsys = 5683 km corresponding to the Heliocentric 
distance oi D — 76 Mpc. In addition, given a estimate for Vsys , the estimate 
for D also depends on the today's Hubble rate value Hq and deceleration pa- 
rameter qo which has been updated [T3| with respect to the values considered 
in these references. The mass modelling of the galaxy is sensitive to several of 
the physical parameters, hence we re-derive these parameters considering the 
converted systemic velocity Vys ~ 5683 km [8l[9] and today's value for the 
Hubble rate Hq = 2.28 x 10^^^ s^^ and deceleration parameter qo = —0.582. 
The derived physical parameters for the galaxy UGC2885 are listed in table [TJ 

Following the modelling setup discussed in [5, to explicitly define the gas 
and stellar matter surface densities on the galaxy plane we are considering the 
neutral Hydrogen Hj 21 cm-line photometric measurements analysed in [7] and 
the K-band photometric measurement analysed in [8]. The gas matter surface 
density is modelled by a thin 2D exponential disk matching the Hj matter 
surface density multiplied by the corrective factor of 1.4 |P and the stellar matter 
surface density is modelled by integrating along the orthogonal spatial direction 
to the galaxy plane the 31? exponential disk matching the K-band luminosity 
in [5]. As for the galaxy bulge we directly fit the excess mass required to the 
rotation curve up to 28 kpc, hence describing the rising of this curve near the 
galactic core. In addition, due to the distinct estimates for D in [7] and [8], it is 
required to scale the surface densities computed in these references accordingly. 
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parameter value ref 



Mtot 


total mass 


1.3 X 10^2 Mq 


i 




mass to light ratio 


4.7Mq/Lq 


s 


MnJMtot 


fractional Hi mass 


« 2.3% 


m 


MoB/Mtot 


disk fractional stellar mass 


« 21% 


m 




rotational velocity at Rmax 


298 (fcms-i) 


mm 


^sys 


measured systemic velocity 


5794 (fcms^i) 


m 




converted systemic velocity 


5683 (fcms^i) 


m 


Rraax 


galaxy semi-major axis 


83.89 (kpc) 


mm 


Z 


cosmological red-shift 


0.01896 


m 


D 


Heliocentric distance 


81.14 Mpc 


m 


H 


Hubble rate 


2.26 X 10-18 


m 


q 


deceleration parameter 


-0.584 


m 



Table 1: Physical parameters of the galaxy UGC2885. The systemic velocity Vsys is 
converted to the heliocentric velocity with respect to the local group and the cosmic 
microwave background 9, and the remaining quantities Rmax, z, D, H and q, are 
computed for the most recent estimates for Ho and qo [13] . 

Hence we obtain 

A'2S 

/in = 0.00791432e~(^7^sfeT^)'5^^j_^jf(ri) (kgm~'^), 

HB = 0.329868 e+T^^^^in^Scnt-osiri) (kgm-^), (^"^^ 

1 ^1 — Rmax 
! Rmax ^ ^1 ^ Rmax ~t- S 
; ^1 ^ Rmax 

where, for a given radial coordinate discretization ri[i] , /is[i] are constant surface 
density disks with edge at ri = (rij^] +ri^i^^)/2, N2S corresponds to the integer 
labelling of the discretized radial coordinate nearest to 28 kpc and the function 
^cut-off describes a smooth cut-off for the galaxy edge near Rmax- 

We now proceed to explicitly define a lattice model for the galaxy UGC2885 
similar to models considered for many body cosmological simulations [23] . Align- 
ing the z axis orthogonally to the galaxy plane such that the galaxy disk lays in 
the plane of constant 9 = tt/2 and considering a regular polar lattice discretiza- 
tion over the plane of the galaxy by imax points along the radial coordinate ri 
and kmax points along the angular coordinate tp we obtain the discrete polar 



3cut-off(?^l) = < 1 ^ 
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coordinates for the lattive points labelled by [i, k] 

~ * ^ J i = 0, 1, . . . , imax j 




i = 
fc = 0,1,. 



■ 5 ^ri 



27r 



(15) 

The lattice faces are cantered at each lattice point [i,k], hence being defined 
as bounded 2D surfaces, and at each lattice point [i, k] we consider a point-like 
mass M[j] of value matching the integrated surface mass-energy density over the 
lattice face[i_fc] 

face 



Mu 



2n 



dririfi{ri) , 



^i[k] 



2 



(16) 



dri ri ^(ri) , 



where fj, is the modelled surface mass density Hence we are considering 

= 1 + imax X kmax point-likc massive bodies over the planar polar lattice. 
Assuming that all the bodies are approximately at a stable circular orbit 
the galaxy model exhibits explicit planar circular symmetry on the plane of 
galaxy such that both the orbital velocities and the gravitational potentials for 
the bodies at each lattice point [i, k] are independent of the angular coordinate 
index k. Hence, noting that the time derivative of the angular coordinate ip 
for each value of the radial coordinate is <Porb[i] = KDrb[i]/''i[i] , the radial 
equations of motion ([T2|) on the galaxy plane are 



fi[j] =0 4=^ Krb[i] = y^-7'iw(-Fjv[i] + Fn^iii) , 



(17) 



and the the potential U and its derivative U' have the following explicit expres- 
sions at each lattice point 

2GMu 



E 



[^] 



[j,fc]#[«,o] 



C/r', = - 



E 

b",fc]/[i,o] 



c^(Ari[,,o]b-,fc])^ 



(18) 



In these expressions we are considering as the test masses the bodies at the 
lattice points [i, 0]. Due to the circular symmetry of the planar model any other 
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value of k for the test masses may generally be considered. As for the functional 
parameter UeS on the galactic plane we consider it to be a radial symmetric func- 
tion such that its values is also independent of the angular lattice index k. In 
the following we will fit this parameter to the experimentally observed galaxy 
rotation curves, hence it is further required to consider a numerical derivative 
to evaluate the derivative of this parameter with respect to the radial coordi- 
nate a'gff. We consider a first order approximation to this derivative. To reduce 
the numerical uncertainty on this approximation we further consider the pa- 
rameterization of the functional parameter by a second order expansion on the 
gravitational field, acs = 3-1- a{l — U)^ . This parameterization does not neces- 
sarily has any physical meaning, it is employed here as a parameterization that 
effectively reduces the relative magnitude of the contribution of the first order 
discrete derivative, hence allowing to numerically fit the functional parameter 
of the model to the existing velocity profiles without requiring to consider an 
higher order approximation to the derivative. Therefore, at each lattice point 
[i, k] on the plane of the galaxy (again independently of the index k) we obtain 
the following expressions for the functional parameter and its derivative with 
respect to the radial coordinate 

aeff[i] 

where we have considered forward finite differences when writing the discrete 
derivative of the functional coefhcient a^jj. This construction allows to describe 
the unknown functional parameter that parameterizes the gravitational inter- 
actions corrections by its values across the discretized lattice points. As for the 
gravitational acceleration contributions to the orbital velocity (|17p are 

Fn[{\ 



^Fh2 [i] 

^rb[i] [13 can now be fitted to the experimental measured rotational velocities. 
At each lattice point the values for the experimentally measured velocity V^xp[i] 
are linearly interpolated for each value of the radial coordinate rj^^j] to the ex- 
perimental profile listed in table 3 of reference [7] . To numerically evaluate the 
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3 + a[i](l - , 



-2t/[',](l -[/[,]) + (l-t/H) 



2 



(19) 



riu+i] - riu] 



1- 



-2{l + q){l-U[l])^-^ 



(20) 



-(l-C/[,])"-w+ilog(l-C/H)a;ff[,] 



values of the functional parameter ags that fit this data we are considering an 
edge smoothing of 5 = 5 kpc and a discretization of the radial coordinate on the 



galaxy plane up to R„ 



.89 kpc with the following lattice properties 



= 74 
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Ar, = 1.20464 (fcpc) 



(21) 



N = 1185 . 

For values of acff[i] ~ —10^ on the galaxy disk the flattening of the galaxy rota- 
tion curve is fully accounted for. The velocity contributions due to the Newton 

y^—ri[i]Fpf[i] and due to the background de- 



gravitational acceleration Vj 



scribed by the El A metric V^ii] = {Fm [i] 

velocity V^rbf, 



A_Ff/2 [j] ) to the total orbital 

^V^j.j + V^j.j are plotted in figure [T] The values for these 

velocities as well as the values of the point-like baryonic masses Mi and the 
coefhcients au] are listed in table [2] in the appendix. We note that reducing 



(a) contributions to the galaxy velocity prolile 



(b) vakies of the I'uiictional parameter (r^rf 




Figure 1: (a) Velocity profiles for UGC2885. The velocity contribution due to the 
Newton acceleration (V/v) is represented by a dashed-dot line, the velocity contribution 
due to ELA metric correction to the gravitational acceleration (Va) is represented by 
a dashed line and the total orbital velocity (Vorb) by a continuous line with the lattice 
points represented by dots; (b) profile for the functional parameter acff on the galactic 
plane, the lattice points are represented by dots. 



the spacing of the lattice, hence increasing the number of the bodies in the 
lattice model and reducing the mass for each of these massive bodies, does not 
change significantly the value of the effective functional parameter aeff- This 
is mainly due to the shift function depending on the total gravitational po- 
tential Un which is approximately maintained constant for each value of the 
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radial coordinate independently of how many discretization points have been 
considered. 

The total baryonic mass for this model is 

Mb = 2.30 X 10" M© . (22) 

However there are contributions to the total galaxy matter due to mass-energy 
density contribution due to the expanding locally anisotropic background correc- 
tions with respect to Schwarzschild backgrounds |16j . In particular it is relevant 
to remark that to ensure that causality is preserved the total mass-energy must 
be strictly positive both on the galaxy disk as well as beyond the galaxy edge. 
For the particular fit discussed here we have considered that beyond the galaxy 
disk edge the value of the parameter acs is approximately a constant ~ 10~^ 
ensuring that the mass-energy density of the background is strictly positive. 
Next we analyse and discuss in detail the mass-energy densities for the galaxy 
model. 



4 Mass-energy density analysis 

To evaluate the mass-energy density pa due to the locally anisotropic back- 
ground within the galaxy plane it is enough to evaluate the expression ((T3)) . 
We have already fitted the functional parameter aoff to the rotation curve on 
the plane of the galaxy as pictured in figure [l] hence it is straight forward to 
evaluate Pa[i,k] at each lattice point [i,k] by excluding the point-like mass at 
this lattice point and consider the contribution of the remaining — 1 point- 
like masses. Noting that the derivatives of the gravitational potential dgll and 
dtpU are null and that pa is circular symmetric on the galaxy plane, the dis- 
crete quantities U[i\, C/j'-j (fTSll . UeS and a'^g (fT9|l already computed are enough 
to actually compute the profile for the mass-energy density Pa[i,o] on the galaxy 
plane. Such profile is plotted in figure [2j 

We can readily verify that there is a negative mass-energy contribution due to 
the expanding locally anisotropic background for radial distances greater than 
21 kpc from the centre of the galaxy. Strictly negative mass-energy densities 
violate causality, therefore are commonly not considered as a physical reality. 
For the specific model discussed here the mass-energy is strictly positive on 
the galaxy plane up to the distance of 89 kpc from the galaxy centre such that 
causality is preserved. To show it explicitly let us note that the surface mass- 
energy contributions pB and po (fT4|) are effective 2D quantities describing the 
3D mass-energy densities ps and pn, respectively. Specifically is fitted to 
a exponential 3D disk in [8] and for po we assume that it is approximately 
described by a Gaussian disk of thickness Ari/2 

PB = 4.60183 X 10^21 g- 11.6367 fcpc- 1.16367 fcpc ^ 

ri / ^ \2 (23) 

Pjj = 2.40691 X 10~^^ e 37.5662 tpc [0.602322 kpc) 

Hence within the galaxy disk the sum of these two contributions plus the con- 
tribution Pa ()13|) is strictly positive. However, beyond the galaxy disk edge 
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ri > 89 kpc, the baryonic mass-energy is approximately null such that the only 
contribution to the total mass-energy density is due to pa- For the galaxy model 
presented here pa > for ri > 89 kpc when UcS > —10^. Also we recall that 
pa does contribute to the total mass of the galaxy in addition to the baryonic 
mass such that to ensure a finite contribution mass for the model we further 
consider a upper radial cut-off Ra above which the functional parameter is null 
aeff = such that explicit spatial isotropy is recovered and the mass-energy 
density above this cut-off exactly matches the one for the expanding Universe 
PRW [H]- Hence both to preserve causality beyond the galaxy edge and to en- 
sure a total finite galaxy mass we consider that the functional parameter acs is 
approximately a constant in-between the radial distance of 89 kpc and Ra 

( -10^ , n e ]88.89 kpc, Ra] , 
aeff = < (24) 

[ , 7-1 e ]i?a,+00[ . 

This ansatz is based only on the assumption that causality must be maintained 
as well as finitude of the total galaxy mass. A more realistic setup and estimate 
for ttoff would require the analysis of direct observational data outside the galaxy 
disk, for instance the gravitational lens effect on the background radiation due 
to the deformation of the gravitational background induced by the galaxy. We 
leave such analysis to another work, here we will proceed our analysis taking 
the ansatz (j24[) as an approximately constant lower cut-off for the functional 
parameter aeff- 

Given the model setup just described, the total mass-energy density is strictly 
positive on the galaxy plane both within the galaxy disk and beyond the galaxy 
edge. The contributions to the total mass-energy density on the galaxy plane 
are plotted in figure [2] 



(a) p„ on the galaxy plane (b) conlributions to p on the galaxy plane 




2(J 40 do SO lOU 120 20 40 60 m lUO 120 

r, (kpc) r, (kpc) 



Figure 2: Mass-Energy profiles for the galaxy model: (a) mass-energy density pa (|13[l 
on the galactic plane, the lattice point are represented by dots; (b) contributions pa 
(thin line) and (dashed line) to the total mass-energy density ptot (thick line) on 
the galaxy plane. The total mass-energy is strictly positive. 

In addition to the analysis of the mass-energy contributions and the ansatz 
for the value of acs beyond the galaxy edge one is led to the question whether 
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our model can further be extended in a 3D neighbourhood of the galaxy, in 
particular whether it is possible to estimate the value of the functional param- 
eter and the mass-energy density along the orthogonal direction to the plan of 
the galaxy. To implement such analysis let us consider the map between met- 
ric (HI) and metric ([7]) and assume the existence of an unique reference profile 
ao corresponding to a reference mass of Mq such that each for the functional 
parameters a„ are defined with respect to ao accordingly to the scaling law ([9]). 
If the 3D profile for ao is known also the 3D profiles for the several a„ 's can be 
computed from the scaling law © and the value for the parameter aeff along 
the orthogonal direction to the galaxy plan can be estimated from the map ([5]) 
by considering some given symmetry for the several functional parameters a„. 

Noting that the galaxy model is explicitly circular symmetric, a discrete 
planar profile ao[i] on the galaxy plane for the functional parameter ao can be 
computed by solving the linear system 



E 



Ml 



0!cH[i] = -1 



[i,fe]#[j,0] 



Mo 



1) log (1 - [/„ 



log l-i7u; 



(25) 



where the potentials Un[i^o]lj,k] correspond to the gravitational potential of the 
point-like massive body located at the lattice point [j, k] evaluated at the lattice 
point [i,0] and Q!o[j,o][i,fc] is the linear interpolated profile for ao[i] evaluated at 



a distance Arj^r, 



from 



"o[o] 



Unli,0]ljM - 



2GM, 



[j] 



»,o]b-,fc] 



Q^O[i,0][i,/c] 



"0[J] 



(Ari[i,o]bMo]-'-i[j])(aom-"o[3+il) 



, (^^i[.,ol[3.fcl-''i[j--il)("o[3~il-"o[j]) 

"o[j-i] + 



nearest(Ari[i^o][j,fc]) 



for Ari[,,o][j,fc] > r^j] 

-''l[j--ll)("0[j^l]-"0[j]) 

A., 

for Ari[i,o][j,fc] < 



(26) 



The potentials and the distances Arij^ o][j,fe] are defined in (fTS]) and j is the 
index of the discrete radial coordinate nearest to Ar^j o][j,fe] ■ 

To ensure that the system of equations (|25l) is not degenerate it is required 
that both the profile aoff[i] and ao[i] have the same length such that the indice 
i for both parameters run from to some Imax- Here we are considering that 
Imax — ^hnax with acff[j>ima^] = ^ 10^. Hcucc, solviug this linear system we 
obtain the profile ao[i] on the galaxy plane. Such profile for our lattice disk 
galaxy model is plotted in figure |3l 

Known the discrete reference profile for ao and some symmetry along the 
orthogonal direction to the galaxy plane it is straight forward to compute the 
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Figure 3: Discrete circular symmetric reference profile ao on tiie galaxy plane. 



value of the functional parameter a^g outside the galaxy plan. Specifically at a 
spatial point orthogonal to the lattice position [i, 0] with a angle 9 with respect 
to the z axis, hence corresponding to the spherical coordinates (ri[j]/ sin(^), 0, 9), 
we obtain 

aeffH(^) = -1 + ^-^^^ ... .... , (27) 

log(l - C/jvliiW) 

where Un[i,o][j,k]{d) is the gravitational potential of the massive body at the 
lattice position [j,k] and UNli]{9) is the total gravitational potential of the A'' 
bodies, both evaluated at the position (ri[j]/ sin(^),O,0) 



2GM, 



b1 



c2 Ari[i,o]b-,fc](6') 



Un[{\{0) = ^Un[ifi]y,k]{6) , (28) 

l3,k] 



Afm^o]U,k]{0) = y(^) -2r-i[i]riy]Cos(¥)[fc])+r2yj . 

In this expressions fi[ifi][j,k]{^) is the Euclidean distance between the lattice 
point [j,k] (i.e. {ri[j],(p[^,9)) and the spatial position {rm^/ sm{9),0,6). 

As for Q:o[i] (9) it is the estimate for the value of the reference parameter ao 
evaluated outside the galaxy plane at the position (ri[i]/ sin(^), 0, ^). To actually 
estimate the value of this parameter at a generic 3D position it is required to 
assume either spherical symmetry of each functional parameter a„ with respect 
to the centre of mass of body n or assuming an explicit anisotropy orthogonally 
to the galaxy plane proportional to (sin^)^'' such that circular symmetry (for 
a fixed value of the coordinate 9) is explicitly maintained. Such anisotropy is 
consistent with the planar model for the galaxy as the matter density distribu- 
tion (a planar disk) explicitly violates spherical symmetry maintaining circular 



19 



symmetry and it allows for a fine-tune of the mass-energy density maintaining 
it strictly positive outside the galaxy plane. In addition we note that as already 
discussed in the introduction, for large radial distances from the galaxy, the 
many body metric must consistently be asymptotically described by one single 
massive body with mass matching the total baryonic mass of the galaxy such 
that spherical symmetry is approximately recovered and the functional param- 
eter ttes is asymptotically independent of both the angular coordinates (p and 
0. Hence we are considering the following assumptions: 

• approximately circular symmetry such that all quantities are independent 
of the angular coordinate (p; 

• anisotropic dependence of the parameter uq on the factor (sin^)^'' which, 
for the particular case of — 0, corresponds to spherical symmetry; 

• for large radial distances the parameter ao is independent of the angular 
coordinate Q approximately matching the value Q!o[/a/ax]- 

Given these assumptions it is straight forward to define the parameter ao eval- 
uated outside the galaxy plane at the position (ri[,j]/ sin(0), 0, 6*) 

- (a\ ^(~- ia\ \ ( ^^l[^o][i■fc](V^) ^^' /oo^ 

aoHlfc*) = "o[/a/ax] + ("o[».o][i,fc](fc') - "o[/max]) — 7^- > (29) 

V ^'^i[i,o][j\fe]lfJ / 

where the sine of the angle 6'[i_o][j,fc] (^) between the orthogonal direction to 
the lattice point [j, fc] and the line between this lattice point and the position 

sin(6'), 0, &) is given by the ratio of the radial distances Afij^ o][j,fc] (7''/2)/Afi[,j o][j,fc] {§) 
and Q:o[i,o][j,fc] (^) is the value of the linearly interpolated reference parameter ao 
corresponding to the distance to the origin of Afi[j o][j,/£] (^)i i-S- the distance to 
ao[o] (EH). 

For = 0, corresponding to spherical symmetric a„'s there exist regions 
outside the galaxy plane where the mass-energy density is negative, while for 
large values of this exponent, > 50000, corresponding to anisotropic a„'s, 
the mass-energy density is strictly positive everywhere. The map for both the 
functional parameter aos and the mass-energy density contribution pa on a 
plane orthogonal to the galaxy containing the galaxy centre is plotted in figurelH 
The several samplings of acfi and pa orthogonal to the galaxy plane up to the 
distance (on the galaxy plane) to the centre of the galaxy of ?'i(gai) = ''i[6i] — 
73A8kpc considered when building the map of figure |4] are plotted in figures [5] 
to [5] in the appendix. 
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a.i) a^ff map for ^^=0 (spherical symmetric a„'s) 



b.i) ffeff map for^e>50000 (anisotropic a^'s) 




a.ii) Pa map for ^£(=0 (spherical symmetric a^'s) 




I 
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''1 (gal) to') 

b.ii) map for ^ti>50000 (anisotropic a„'s) 



<-l,7 X lU"^^ ■ 

-100 -50 50 100 
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Figure 4: The map of the functional parameter cceff and the mass-energy density pa 
on a plane orthogonal to the galaxy containing the center of the galaxy, the x axis 
corresponds to the galactic plane such that ri(gai) is the distance on the galaxy plane 
to the centre of the galaxy and the y axis corresponds to the orthogonal distance to 
the plane of the galaxy: a) assuming spherical symmetry, = 0, for all q„'s there 
exist regions outside the galaxy plane where the mass-energy is negative; b) assuming 
anisotropy along the orthogonal direction to the galaxy plane, for > 50000, the 
mass-energy density is strictly positive outside the galaxy plane. 
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Noting that, independently of the direction, for large radial distances to 
the centre of the galaxy the functional parameter is approximately constant 
aeff = —10^ up to the upper cut-off Ra (1^ . it is straight forward to estimate 
the relative mass-energy contribution for the relative cosmological mass-energy 
density corresponding to the galaxy model developed in this work. Let us con- 
sider a relative large distance from the galaxy such that the gravitational field of 
the galaxy is approximately given by one single point-like massive body, hence 
the mass contribution in addition to the expanding Universe background mass- 
energy density pRw is 

M„ = 47r r\{p^ - pj,w)dr^ ^^V^' UiRcT^"^""-^ ' ^) RI ■ (30) 

Assuming that all the baryonic matter has a similar contribution to the relative 
cosmological mass-energy contribution Ma ~ 4.96 Mf, [13] we obtain for our 
galaxy model (f22|) 

Mtot = 13.68 X lO^A^Q , M„ = 11.38 x lO" Mq ^ i?„ = 9574.61 kpc (31) 

Recalling that within the ACDM cosmological model such contribution must 
be included in the relative mass-energy density of dark matter, the extended 
background described by the ELA metric up to Ra is interpreted as the heuristic 
dark matter halo commonly employed to describe dark matter effects within 
galaxies. 

5 Discussion and Outlook 

By fitting the expanding locally anisotropic metric functional parameter acfi ([4]) 
we have fully describe the fiattening of the rotation curve for the galaxy UGC2885 
Based on a lattice model we have estimate values of this functional parameter 
of order ~ —10"'' within the galaxy disk simultaneously maintaining compati- 
bility with the well established short-scale gravitational laws (the Schwarzschild 
metric is approximately recovered at planetary scales) and cosmological scale 
gravitational laws (the Universe expansion). Phenomenologically this is inter- 
preted as a unaccounted gravitational interaction acting on baryonic matter that 
reproduces the effects attributed to cold dark matter, namely the flattening of 
galaxies rotation curves. 

Also we have analysed the mass-energy density on the galaxy plane showing 
that although there are negative contributions from the extended background 
described by the ELA metric this quantity is strictly positive within the galaxy 
disk. Outside the galaxy disk, both on the plane of the galaxy and outside the 
galaxy plane, to ensure that the mass-energy density is strictly positive, it is 
considered a lower absolute value for the functional parameter (specifically ~ 
— 10~^). This construction implies that the metric functional parameter is larger 
in absolute value in the vicinity of baryonic matter, i.e. on the galaxy plane, 
than in the vicinity of the galaxy such that the gravitational corrections are 
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localized in the galaxy disk. In addition the functional parameter is null above 
the upper cut-off Ra = 9574.61 kpc such that the ratio of the galaxy baryonic 
matter Mb to Ma , the mass contribution of the extended background in addition 
to the cosmological background matches the ratio of the baryonic to cold dark 
matter relative cosmological mass-energy densities Jlb/flc = Mb/M^ = 4.96 
consistently with the physical interpretation that the ELA metric constitutes a 
parameterization of cold dark matter. 

A relevant question raised by the construction developed here is whether the 
gravitational corrections due to the expanding locally anisotropic background 
to the redshift of radiation emitted from the galaxy are significant or not [TB]. 
To answer this question let us recall that when computing the profile for the 
functional parameter we have excluded the point-like masses at each lattice 
point as these are considered test masses on the background of the remaining 
— 1 massive bodies on the galaxy. When considering an external point mass 
in the vicinity of each body on the galaxy, the effective functional parameter 
for the A^ body system, as perceived by the external test mass, must decrease 
in absolute value as near massive objects the Schwarzschild metric is asymp- 
totically recover. Hence for a external test mass in the galaxy plane, near each 
massive body the functional parameter decreases in absolute value such that the 
red-shift corrections are negligible for radiation emitted from each body on the 
galaxy. As an analogy let us note that for many-body gravitational systems on 
Schwarzschild geometries, when computing orbits a given planet is considered a 
test mass such that its own gravitational field is not accounted in the equations 
of motion, however the frequency-shift for radiation emitted from each planet 
is mainly accounted for by the planet gravitational field. 

As future directions of research we note that to develop a more detailed 
model for galaxies parameterized by a many body expanding locally anisotropic 
background it must be derived a continuous description of the intergalactic 
matter (gas) not contained on stars. Also, within a more detailed model, the 
assumption that all masses have a similar profile for the functional parameter 
given by the scaling law © may be an over simplification as the localization of 
the gravitational corrections in the vicinity of galactic plane indicates that the 
functional parameter for each of the many bodies on the galaxy may depend 
on the remaining masses in the gravitational system, hence not being exactly 
circular symmetric as considered in the simplified model developed in this work. 

As a final remark let us note that a direct analysis in the neighbourhood of 
galaxies which may either confirm or dismiss the construction suggested here is 
achievable, for instance, by analysing the lens effect of the background radiation 
near galaxies. Also within the Solar System a test mass orbiting orthogonally 
to the planetary system plane would allow for a direct measurement of the 
localization of the gravitational corrections encoded in the expanding locally 
anisotropic metric. We leave these analysis to another work. 
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A Appendix 

In this appendix we list the values for the galaxy model masses, functional 
parameter and velocity contributions corresponding to figure[T]in tableland the 
samplings along the orthogonal to the galaxy plane of the functional parameter 
and mass-energy densities considered to plot the maps of figure S] in figures [S] 
toU 
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Table 2: Values of the point-like masses Mjj], functional parameter QefF[i], and velocity 
contributions yjv[i] and Va[i] to Vorb[i]- 
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a„(z)torr,„=1.204M 



a,„(z)forr,,ja,=4.S1858 



r,„(z) for r, ^,=8.43251 



a,i,(z)torr,|ga,=12.04M 



actf(z) for ri(gai)=15.6604 



i>dt(z)fornw=19.2743 




0M2) for n ,gai)=22.8882 




ieff(z)forn(gai)=26.5022 



i,„(z)forr,^l=30.1161 




Zefr(z}forn(gai)=33.73 



«,„(z)forn„=37.344 



i,tt(z) for n (^,=40.9579 






ff,„(z)forni,a,=44.571S 




t,ii(z)forn(gai)=4S.1858 



i.i,(z)forr,„=51.7997 




i,(z) for r, |,a,=55.4136 




.„(z)forr|,ga,=59.0276 



r,„(z) for r, ^,=62.6415 



lerf(z) for n ,gai)=66.2554 



.,i(z)forn(gai)=69,S694 



i.i,(z)forr,„=73.4833 



Figure 5: Estimates for Qefi outside of tlie galaxy plane assuming spherical symmetry 
of the functional parameters a„'s {^e = 0). 



27 



(>„(z)torriju,=1.20464 



;::a 



p„(z)forri|gai,=4.81858 



p„(z)for r,,gj)=8.43251 



p„(z)forri(gai)=12.0464 



Pa{-£) for r\ (gai)=15.6604 



Pa(z)forri(gii)=19.2743 



Ptt(z)forri(gai)=22.8882 



/3a{z) for ri ,gai)=26.5022 



Pa(z) for ri(gai)=30.1161 





p„(z) for n (gai)=33.73 



/3a(z) for ri |gai)=37.344 



Pa(z)forri(ga,)=40.9579 



Pa(z) for n (gai)=44.5718 



/3tt{z)forri(£ai)=48.1858 



Pa(z)forri(gii)=51,7997 






p„(z)forrng^l=55.4136 



p„(z) for /-i ^,=59.0276 



p„(z)forri,gj)=62,6415 



Ptt(z) for ri (gall =66.2554 




Pa(i) for ri ,gai)=69.8694 



/[3„(z)fOTri(gii)=73.4833 



Figure 6: Estimates for outside of the galaxy plane assuming spherical symmetry 
of the functional parameters a^'s (^^ = 0). 
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a^(z) tor r, (^1=1.20464 



(iai(z)forr,„=4.81S5S 



a,,(z)torr,„=S,43251 




200 250 



aeff(z) for ri (gai|=12.0464 



OTeff (z) for ri (gai)=15.6604 



* -0., 



<»„(z) for r, (^1=19.2743 



ffeff (z) for ri |gai)=22.8i 



aett(z) for ri (gai)=26.5022 



iz«,(z) for r, (5,1=30.1161 



«eff (z) for ri (gai)=33.73 



ff,j,(z)forr,(sa,=37.344 



tteff{z)forfiM=40.9579 



o«,(z)fori-,ua,=44.571S 



ae£f(z)forri,gai)=4S.1858 



50 100 




a.j(z)forr,(gj|=51.7997 



o„(z)forr,„=55.4136 



««i(z)forr,„=59.0276 



ffeifCz) for fi ,gai,=62.6415 



aeff (z) for ri |gai)=66.2554 



aeff(z)fOTri(gai)=69,Si 



Odi(z)forr,(,j|=73.4833 



Figure 7: Estimates for aes outside of tlie galaxy plane assuming anisotropy of the 
functional parameters a„'s {^e > 50000). 
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(>,(z)forr,„,=1.20464 



/>,(z)torr,j|j,=4.81858 



/j„(z)tori-m)=8.43251 




